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Strong coupling expansion of chiral models 
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A general precedure is outlined for an algorithmic implementation of the strong coupling expansion of lattice 
chiral models on arbitrary lattices. A symbolic character expansion in terms of connected values of group integrals 
on skeleton diagrams may be obtained by a fully computerized approach. 



1. Introduction. 

The main motivations for a strong-coupling ap- 
proach to asymptotically free matrix-valued field 
theories (like 2-d principal chiral models) are the 
following: 

(i) Search for a scaling region, within the con- 
vergence radius of the strong coupling series, 
where the continuum results may be extracted 
by appropriate resummation. Numerical evidence 
favors this possibility [1]. 

(ii) Extension of the scaling region, essentially 
by adopting the "energy" scheme [3] in the defini- 
tion of the coupling directly in the strong coupling 
series and studying convergence and scaling pro- 
perties. Numerical evidence indicates asymptotic 
scaling within the strong coupling domain [1,2]. 

(iii) Study of the relevance of the planar li- 
mit N — > oo, which in turn leads to an analysis 
of the double-scaling limit properties of many- 
matrix systems, in particular a determination of 
the location and nature of the criticality and eva- 
luation of the relevant critical exponents. Nume- 
rical evidence points at a second order transition 
at /3 ~ 0.306 for the 2-d chiral models on the 
square lattice [2]. 

2. Strong coupling expansion. 

The strong coupling expansion is best perfor- 
med by adopting the character expansion of the 
action [4,5] 

e 2N,K^U = e *'*(«£ d(r)Z(r)X(r)(E0> (1) 

(r) 
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where F(/3) is the free-energy of the single-link 
model, X(r) are the characters and d( r ) the dimen- 
sions of the irreducible representations of U(N) 
and Z( r ) are the character coefficients and are 
naturally organized into a strong-coupling series. 
Notable simplifications occur in the large- TV limit 
of the strong coupling expansion: we must howe- 
ver notice that strong coupling and large- TV limit 
are not strictly commuting, which implies that 
some care is needed in the extrapolations, that 
are however supported by numerical evidence. 

In order to sketch the logic of an algorithmic 
implementation of the strong coupling expansion 
we must introduce a few definitions. Lattice 
functional integrals are reduced by the character 
expansion to sums of group integrals on nontrivial 
assignements of representations on the lattice li- 
nks. Assignments can be grouped into "oriented 
configurations" characterized only by the number 
of "quarks" n + and "antiquarks" n_ running on 
each lattice link. Configurations (n + + n_ = n 
constant) can be turned into skeleton diagrams, 
where each (abstract) link is characterized by 
n and the length of the associated lattice path. 
Configurations sharing the same skeleton have the 
same value of the group integral. We therefore 
define the "geometrical factor" as the number of 
configurations on each skeleton, and the "group- 
theoretical factor" as the value of the group inte- 
gral for each skeleton. 

The geometrical factor can be evaluated by a 
fully computerized procedure: 

- generate all nonbacktracking random walks 
with fixed length, i.e. fixed order in strong cou- 
pling; 
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- keep one and only one walk for each different 
configuration; 

- identify the skeletons of configurations and 
compare them to each other implementing per- 
mutation symmetries; 

- collect all configurations with equivalent ske- 
leton. 

The group-theoretical factor requires compu- 
ting group integrals, and is in principle a com- 
pletely solved problem, but in practice it re- 
quires generating algorithmically the Clebsch- 
Gordan coefficients for the decomposition of pro- 
ducts of arbitrary irreducible representations of 
U(N). Some progress in the evaluation of group- 
theoretical factors is obtained by eliminating bub- 
bles in skeletons diagrams and reducing the pro- 
blem to evaluating only superskeletons group in- 
tegrals (a much smaller number of objects): su- 
perskeletons are diagrams where no multiple con- 
nections between sites are allowed; they are ob- 
tained by use of the simpler decomposition of pro- 
ducts of characters. 

By known results of group integration [6], all in- 
tegrals relevant to rather high orders can be expli- 
citly computed. By exploiting geometric proper- 
ties of oriented configurations, we may evaluate 
the connected contribution of each skeleton to the 
expectation value of a physical quantity. Connec- 
ted group-theoretical factors are called "poten- 
tials": thay are the quantities admitting a finite 
large- TV limit. Wide families of potentials, rela- 
ted to the same superskeleton, can be computed 
symbolically in terms of the Z( r ) and numerically 
in terms of N and /3. The symbolic character 
expansion in terms of potentials is obtained by a 
fully computerized process. 

3. The fundamental Green's function. 

A few important considerations concern 
Green's functions 

G(x) = l<Tr[t^M0)t]>, 

G(p) = X>**G(z). (2) 

X 

Truncated strong coupling expansions admitting 
a reinterpretation as summations over paths can 



be see as originated by generalized Gaussian mo- 
dels with non-nearest interactions. At q th order 
in strong coupling expansion the inverse propa- 
gator admit the expansion 

Gip)- 1 = A (p) + A^p 2 

+E E A uA P)P 2s [( P 2 ) 2 -pf~ s %) 

u=2 s =0 

u — s=even 

where the number of independent coefficients 
A u , s ~ ft Zu is equal to the number of indepen- 
dent effective couplings. Strong coupling series 
for A u>s are an efficient way of storing all the 
strong coupling information concerning G{x). 

Physical quantities may be extracted directly 
from A u<s : 

- magnetic susceptibility: x = l/^o- 

- second-moment correlation length: {x 2 )q = 
Ai/Ao. 

- wavefunction renormalization: Zq = 1/Ai. 

- True mass gap jj, s by inversion of 

x = 2(1 - eh//,) , (4) 
where x is the perturbative solution of 

oo 

A + A lX + ^ A UtU x u = 0. (5) 

u = 2 

Thermodynamical quantities like energy and 
specific heat may also easily be extracted. 
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